Synchrotron-source micro-x-ray computed tomography for examining butterfly eyes.

Comparative anatomy is an important tool for investigating evolutionary relationships among species, but the lack of scalable imaging tools and stains for rapidly mapping the microscale anatomies of related species poses a major impediment to using comparative anatomy approaches for identifying evolutionary adaptations. We describe a method using synchrotron source micro-x-ray computed tomography (syn-µXCT) combined with machine learning algorithms for high-throughput imaging of Lepidoptera (i.e., butterfly and moth) eyes. Our pipeline allows for imaging at rates of ~15 min/mm3 at 600 nm3 resolution. Image contrast is generated using standard electron microscopy labeling approaches (e.g., osmium tetroxide) that unbiasedly labels all cellular membranes in a species-independent manner thus removing any barrier to imaging any species of interest. To demonstrate the power of the method, we analyzed the 3D morphologies of butterfly crystalline cones, a part of the visual system associated with acuity and sensitivity and found significant variation within six butterfly individuals. Despite this variation, a classic measure of optimization, the ratio of interommatidial angle to resolving power of ommatidia, largely agrees with early work on eye geometry across species. We show that this method can successfully be used to determine compound eye organization and crystalline cone morphology. Our novel pipeline provides for fast, scalable visualization and analysis of eye anatomies that can be applied to any arthropod species, enabling new questions about evolutionary adaptations of compound eyes and beyond.

Supergene evolution via gain of auto-regulation.

The development of complex phenotypes requires the coordinated action of many genes across space and time, yet many species have evolved the ability to develop multiple discrete, alternate phenotypes1-3. Such polymorphisms are often controlled by supergenes, sets of tightly-linked mutations in one or more loci that function together to produce a complex phenotype4. Although theories of supergene evolution are well-established, the mutations that cause functional differences between supergene alleles remain essentially unknown. doublesex is the master regulator of insect sexual differentiation but functions as a supergene in multiple Papilio swallowtail butterflies, where divergent dsx alleles control development of discrete non-mimetic or mimetic female wing color patterns5-7. Here we demonstrate that the functional elements of the mimetic allele in Papilio alphenor are six new cis-regulatory elements (CREs) spread across 150 kb that are bound by DSX itself. Our findings provide experimental support to classic supergene theory and suggest that the evolution of auto-regulation may provide a simple route to supergene origination and to the co-option of pleiotropic genes into new developmental roles.

Acute and Long-Term Consequences of Co-opted doublesex on the Development of Mimetic Butterfly Color Patterns.

Novel phenotypes are increasingly recognized to have evolved by co-option of conserved genes into new developmental contexts, yet the process by which co-opted genes modify existing developmental programs remains obscure. Here, we provide insight into this process by characterizing the role of co-opted doublesex in butterfly wing color pattern development. dsx is the master regulator of insect sex differentiation but has been co-opted to control the switch between discrete nonmimetic and mimetic patterns in Papilio alphenor and its relatives through the evolution of novel mimetic alleles. We found dynamic spatial and temporal expression pattern differences between mimetic and nonmimetic butterflies throughout wing development. A mimetic color pattern program is switched on by a pulse of dsx expression in early pupal development that causes acute and long-term differential gene expression, particularly in Wnt and Hedgehog signaling pathways. RNAi suggested opposing, novel roles for these pathways in mimetic pattern development. Importantly, Dsx co-option caused Engrailed, a primary target of Hedgehog signaling, to gain a novel expression domain early in pupal wing development that is propagated through mid-pupal development to specify novel mimetic patterns despite becoming decoupled from Dsx expression itself. Altogether, our findings provide multiple views into how co-opted genes can both cause and elicit changes to conserved networks and pathways to result in development of novel, adaptive phenotypes.

New molecular insights into butterfly pigmentation.

Hanly et al.1 and Nishida et al.2 use distinct approaches to provide exceptional lessons regarding the genetic, molecular, morphological, and biochemical bases of butterfly wing pigmentation. These mechanistic insights collectively have important implications for our understanding of phenotype evolution.

Migration genetics take flight: genetic and genomic insights into monarch butterfly migration.

Monarch butterflies have emerged as a model system in migration genetics. Despite inherent challenges associated with studying the integrative phenotypes that characterize migration, recent research has highlighted genes and transcriptional networks underlying aspects of the monarch's migratory syndrome. Circadian clock genes and the vitamin A synthesis pathway regulate reproductive diapause initiation, while diapause termination appears to involve calcium and insulin signaling. Comparative approaches have highlighted genes that distinguish migratory and nonmigratory monarch populations, as well as genes associated with natural variation in propensity to initiate diapause. Population genetic techniques demonstrate that seasonal migration can collapse patterns of spatial structure at continental scales, whereas loss of migration can drive differentiation between even nearby populations. Finally, population genetics can be applied to reconstruct the monarch's evolutionary history and search for contemporary demographic changes, which can provide relevant context for understanding recently observed declines in overwintering North American monarch numbers.

Divergent expression of aristaless1 and aristaless2 during embryonic appendage and pupal wing development in butterflies.

BACKGROUND: Gene duplication events are critical for the evolution of new gene functions. Aristaless is a major regulator of distinct developmental processes. It is most known for its role during appendage development across animals. However, more recently other distinct biological functions have been described for this gene and its duplicates. Butterflies and moths have two copies of aristaless, aristaless1 (al1) and aristaless2 (al2), as a result of a gene duplication event. Previous work in Heliconius has shown that both copies appear to have novel functions related to wing color patterning. Here we expand our knowledge of the expression profiles associated with both ancestral and novel functions of Al1 across embryogenesis and wing pigmentation. Furthermore, we characterize Al2 expression, providing a comparative framework between gene copies within the same species, allowing us to understand the origin of new functions following gene duplication. RESULTS: Our work shows that the expression of both Al1 and Al2 is associated with the ancestral function of sensory appendage (leg, mouth, spines, and eyes) development in embryos. Interestingly, Al1 exhibits higher expression earlier in embryogenesis while the highest levels of Al2 expression are shifted to later stages of embryonic development. Furthermore, Al1 localization appears extranuclear while Al2 co-localizes tightly with nuclei earlier, and then also expands outside the nucleus later in development. Cellular expression of Al1 and Al2 in pupal wings is broadly consistent with patterns observed during embryogenesis. We also describe, for the first time, how Al1 localization appears to correlate with zones of anterior/posterior elongation of the body during embryonic growth, showcasing a possible new function related to Aristaless' previously described role in appendage extension. CONCLUSIONS: Overall, our data suggest that while both gene copies play a role in embryogenesis and wing pigmentation, the duplicates have diverged temporally and mechanistically across those functions. Our study helps clarify principles behind sub-functionalization and gene expression evolution associated with developmental functions following gene duplication events.

aristaless1 has a dual role in appendage formation and wing color specification during butterfly development.

BACKGROUND: Highly diverse butterfly wing patterns have emerged as a powerful system for understanding the genetic basis of phenotypic variation. While the genetic basis of this pattern variation is being clarified, the precise developmental pathways linking genotype to phenotype are not well understood. The gene aristaless, which plays a role in appendage patterning and extension, has been duplicated in Lepidoptera. One copy, aristaless1, has been shown to control a white/yellow color switch in the butterfly Heliconius cydno, suggesting a novel function associated with color patterning and pigmentation. Here we investigate the developmental basis of al1 in embryos, larvae, and pupae using new antibodies, CRISPR/Cas9, RNAi, qPCR assays of downstream targets, and pharmacological manipulation of an upstream activator. RESULTS: We find that Al1 is expressed at the distal tips of developing embryonic appendages consistent with its ancestral role. In developing wings, we observe Al1 accumulation within developing scale cells of white H. cydno during early pupation while yellow scale cells exhibit little Al1 at this time point. Reduced Al1 expression is also associated with yellow scale development in al1 knockouts and knockdowns. We propose that Al1 expression in future white scales might be related to an observed downregulation of the enzyme Cinnabar and other genes that synthesize and transport the yellow pigment, 3-hydroxykynurenine (3-OHK). Finally, we provide evidence that Al1 activation is under the control of Wnt signaling. CONCLUSIONS: We propose a model in which high levels of Al1 during early pupation, which are mediated by Wnt, are important for melanic pigmentation and specifying white portions of the wing while reduced levels of Al1 during early pupation promote upregulation of proteins needed to move and synthesize 3-OHK, promoting yellow pigmentation. In addition, we discuss how the ancestral role of aristaless in appendage extension may be relevant in understanding the cellular mechanism behind color patterning in the context of the heterochrony hypothesis.

Seasonal plasticity in morphology and metabolism differs between migratory North American and resident Costa Rican monarch butterflies.

Environmental heterogeneity in temperate latitudes is expected to maintain seasonally plastic life-history strategies that include the tuning of morphologies and metabolism that support overwintering. For species that have expanded their ranges into tropical latitudes, it is unclear the extent to which the capacity for plasticity will be maintained or will erode with disuse. The migratory generations of the North American (NA) monarch butterfly Danaus plexippus lead distinctly different lives from their summer generation NA parents and their tropical descendants living in Costa Rica (CR). NA migratory monarchs postpone reproduction, travel thousands of kilometers south to overwinter in Mexico, and subsist on little food for months. Whether recently dispersed populations of monarchs such as those in Costa Rica, which are no longer subject to selection imposed by migration, retain ancestral seasonal plasticity is unclear. To investigate the differences in seasonal plasticity, we reared the NA and CR monarchs in summer and autumn in Illinois, USA, and measured the seasonal reaction norms for aspects of morphology and metabolism related to flight. NA monarchs were seasonally plastic in forewing and thorax size, increasing wing area and thorax to body mass ratio in autumn. While CR monarchs increased thorax mass in autumn, they did not increase the area of the forewing. NA monarchs maintained similar resting and maximal flight metabolic rates across seasons. However, CR monarchs had elevated metabolic rates in autumn. Our findings suggest that the recent expansion of monarchs into habitats that support year-round breeding may be accompanied by (1) the loss of some aspects of morphological plasticity as well as (2) the underlying physiological mechanisms that maintain metabolic homeostasis in the face of temperature heterogeneity.

Museum genomics reveals the Xerces blue butterfly (Glaucopsyche xerces) was a distinct species driven to extinction.

The last Xerces blue butterfly was seen in the early 1940s, and its extinction is credited to human urban development. This butterfly has become a North American icon for insect conservation, but some have questioned whether it was truly a distinct species, or simply an isolated population of another living species. To address this question, we leveraged next-generation sequencing using a 93-year-old museum specimen. We applied a genome skimming strategy that aimed for the organellar genome and high-copy fractions of the nuclear genome by a shallow sequencing approach. From these data, we were able to recover over 200 million nucleotides, which assembled into several phylogenetically informative markers and the near-complete mitochondrial genome. From our phylogenetic analyses and haplotype network analysis we conclude that the Xerces blue butterfly was a distinct species driven to extinction.

The evolution and genetics of sexually dimorphic 'dual' mimicry in the butterfly Elymnias hypermnestra.

Sexual dimorphism is a major component of morphological variation across the tree of life, but the mechanisms underlying phenotypic differences between sexes of a single species are poorly understood. We examined the population genomics and biogeography of the common palmfly Elymnias hypermnestra, a dual mimic in which female wing colour patterns are either dark brown (melanic) or bright orange, mimicking toxic Euploea and Danaus species, respectively. As males always have a melanic wing colour pattern, this makes E. hypermnestra a fascinating model organism in which populations vary in sexual dimorphism. Population structure analysis revealed that there were three genetically distinct E. hypermnestra populations, which we further validated by creating a phylogenomic species tree and inferring historical barriers to gene flow. This species tree demonstrated that multiple lineages with orange females do not form a monophyletic group, and the same is true of clades with melanic females. We identified two single nucleotide polymorphisms (SNPs) near the colour patterning gene WntA that were significantly associated with the female colour pattern polymorphism, suggesting that this gene affects sexual dimorphism. Given WntA's role in colour patterning across Nymphalidae, E. hypermnestra females demonstrate the repeatability of the evolution of sexual dimorphism.

Migration behaviour of commercial monarchs reared outdoors and wild-derived monarchs reared indoors.

Captive rearing of monarch butterflies is a commercial and personal pursuit enjoyed by many different groups and individuals. However, the practice remains controversial, especially after new evidence showed that both a group of commercially derived monarchs reared outdoors and a group of wild-derived but indoor-reared monarchs failed to orient south, unlike wild-derived monarchs reared outdoors. To more fully characterize the mechanisms responsible for the loss of orientation in both commercial and indoor-reared monarchs, we performed flight simulator experiments to determine (i) whether any fraction of commercial monarchs maintains a southern heading over multiple tests, and (ii) whether indoor conditions with the addition of sunlight can induce southern flight in wild-derived monarchs. Commercial monarchs changed their flight direction more often over the course of multiple tests than wild-derived monarchs. While as a group the commercial monarchs did not fly south on average, a subset of individuals did orient south over multiple tests, potentially explaining the discordance between flight simulator assays and the recovery of tagged commercial monarchs at overwintering locations. We also show that even when raised indoors with sunlight, wild-derived monarchs did not consistently orient south in the flight simulator, though wild-derived monarchs reared outdoors did orient south.

Divergence, gene flow, and the origin of leapfrog geographic distributions: The history of colour pattern variation in Phyllobates poison-dart frogs.

The geographic distribution of phenotypic variation among closely related populations is a valuable source of information about the evolutionary processes that generate and maintain biodiversity. Leapfrog distributions, in which phenotypically similar populations are disjunctly distributed and separated by one or more phenotypically distinct populations, represent geographic replicates for the existence of a phenotype, and are therefore especially informative. These geographic patterns have mostly been studied from phylogenetic perspectives to understand how common ancestry and divergent evolution drive their formation. Other processes, such as gene flow between populations, have not received as much attention. Here, we investigate the roles of divergence and gene flow between populations in the origin and maintenance of a leapfrog distribution in Phyllobates poison frogs. We found evidence for high levels of gene flow between neighbouring populations but not over long distances, indicating that gene flow between populations exhibiting the central phenotype may have a homogenizing effect that maintains their similarity, and that introgression between 'leapfroging' taxa has not played a prominent role as a driver of phenotypic diversity in Phyllobates. Although phylogenetic analyses suggest that the leapfrog distribution was formed through independent evolution of the peripheral (i.e. leapfrogging) populations, the elevated levels of gene flow between geographically close populations poise alternative scenarios, such as the history of phenotypic change becoming decoupled from genome-averaged patterns of divergence, which we cannot rule out. These results highlight the importance of incorporating gene flow between populations into the study of geographic variation in phenotypes, both as a driver of phenotypic diversity and as a confounding factor of phylogeographic inferences.

The roles of hybridization and habitat fragmentation in the evolution of Brazil's enigmatic longwing butterflies, Heliconius nattereri and H. hermathena.

BACKGROUND: Heliconius butterflies are widely distributed across the Neotropics and have evolved a stunning array of wing color patterns that mediate Müllerian mimicry and mating behavior. Their rapid radiation has been strongly influenced by hybridization, which has created new species and allowed sharing of color patterning alleles between mimetic species pairs. While these processes have frequently been observed in widespread species with contiguous distributions, many Heliconius species inhabit patchy or rare habitats that may strongly influence the origin and spread of species and color patterns. Here, we assess the effects of historical population fragmentation and unique biology on the origins, genetic health, and color pattern evolution of two rare and sparsely distributed Brazilian butterflies, Heliconius hermathena and Heliconius nattereri. RESULTS: We assembled genomes and re-sequenced whole genomes of eight H. nattereri and 71 H. hermathena individuals. These species harbor little genetic diversity, skewed site frequency spectra, and high deleterious mutation loads consistent with recent population bottlenecks. Heliconius hermathena consists of discrete, strongly isolated populations that likely arose from a single population that dispersed after the last glacial maximum. Despite having a unique color pattern combination that suggested a hybrid origin, we found no genome-wide evidence that H. hermathena is a hybrid species. However, H. hermathena mimicry evolved via introgression, from co-mimetic Heliconius erato, of a small genomic region upstream of the color patterning gene cortex. CONCLUSIONS: Heliconius hermathena and H. nattereri population fragmentation, potentially driven by historical climate change and recent deforestation, has significantly reduced the genetic health of these rare species. Our results contribute to a growing body of evidence that introgression of color patterning alleles between co-mimetic species appears to be a general feature of Heliconius evolution.

Genomic evidence for gene flow between monarchs with divergent migratory phenotypes and flight performance.

Monarch butterflies are known for their spectacular annual migration in eastern North America, with millions of monarchs flying up to 4,500 km to overwintering sites in central Mexico. Monarchs also live west of the Rocky Mountains, where they travel shorter distances to overwinter along the Pacific Coast. It is often assumed that eastern and western monarchs form distinct evolutionary units, but genomic studies to support this notion are lacking. We used a tethered flight mill to show that migratory eastern monarchs have greater flight performance than western monarchs, consistent with their greater migratory distances. However, analysing more than 20 million SNPs in 43 monarch genomes, we found no evidence for genomic differentiation between eastern and western monarchs. Genomic analysis also showed identical and low levels of genetic diversity, and demographic analyses indicated similar effective population sizes and ongoing gene flow between eastern and western monarchs. Gene expression analysis of a subset of candidate genes during active flight revealed differential gene expression related to nonmuscular motor activity. Our results demonstrate that eastern and western monarchs maintain migratory differences despite ongoing gene flow, and suggest that migratory differences between eastern and western monarchs are not driven by select major-effects alleles. Instead, variation in migratory distance and destination may be driven by environmentally induced differential gene expression or by many alleles of small effect.

Disentangling Population History and Character Evolution among Hybridizing Lineages.

Understanding the origin and maintenance of adaptive phenotypic novelty is a central goal of evolutionary biology. However, both hybridization and incomplete lineage sorting can lead to genealogical discordance between the regions of the genome underlying adaptive traits and the remainder of the genome, decoupling inferences about character evolution from population history. Here, to disentangle these effects, we investigated the evolutionary origins and maintenance of Batesian mimicry between North American admiral butterflies (Limenitis arthemis) and their chemically defended model (Battus philenor) using a combination of de novo genome sequencing, whole-genome resequencing, and statistical introgression mapping. Our results suggest that balancing selection, arising from geographic variation in the presence or absence of the unpalatable model, has maintained two deeply divergent color patterning haplotypes that have been repeatedly sieved among distinct mimetic and nonmimetic lineages of Limenitis via introgressive hybridization.

A shared genetic basis of mimicry across swallowtail butterflies points to ancestral co-option of doublesex.

Uncovering whether convergent adaptations share a genetic basis is consequential for understanding the evolution of phenotypic diversity. This information can help us understand the extent to which shared ancestry or independent evolution shape adaptive phenotypes. In this study, we first ask whether the same genes underlie polymorphic mimicry in Papilio swallowtail butterflies. By comparing signatures of genetic variation between polymorphic and monomorphic species, we then investigate how ancestral variation, hybridization, and independent evolution contributed to wing pattern diversity in this group. We report that a single gene, doublesex (dsx), controls mimicry across multiple taxa, but with species-specific patterns of genetic differentiation and linkage disequilibrium. In contrast to widespread examples of phenotypic evolution driven by introgression, our analyses reveal distinct mimicry alleles. We conclude that mimicry evolution in this group was likely facilitated by ancestral polymorphism resulting from early co-option of dsx as a mimicry locus, and that evolutionary turnover of dsx alleles may underlie the wing pattern diversity of extant polymorphic and monomorphic lineages.

Genomic architecture and introgression shape a butterfly radiation.

We used 20 de novo genome assemblies to probe the speciation history and architecture of gene flow in rapidly radiating Heliconius butterflies. Our tests to distinguish incomplete lineage sorting from introgression indicate that gene flow has obscured several ancient phylogenetic relationships in this group over large swathes of the genome. Introgressed loci are underrepresented in low-recombination and gene-rich regions, consistent with the purging of foreign alleles more tightly linked to incompatibility loci. Here, we identify a hitherto unknown inversion that traps a color pattern switch locus. We infer that this inversion was transferred between lineages by introgression and is convergent with a similar rearrangement in another part of the genus. These multiple de novo genome sequences enable improved understanding of the importance of introgression and selective processes in adaptive radiation.

Comparative Transcriptomics Provides Insights into Reticulate and Adaptive Evolution of a Butterfly Radiation.

Butterfly eyes are complex organs that are composed of a diversity of proteins and they play a central role in visual signaling and ultimately, speciation, and adaptation. Here, we utilized the whole eye transcriptome to obtain a more holistic view of the evolution of the butterfly eye while accounting for speciation events that co-occur with ancient hybridization. We sequenced and assembled transcriptomes from adult female eyes of eight species representing all major clades of the Heliconius genus and an additional outgroup species, Dryas iulia. We identified 4,042 orthologous genes shared across all transcriptome data sets and constructed a transcriptome-wide phylogeny, which revealed topological discordance with the mitochondrial phylogenetic tree in the Heliconius pupal mating clade. We then estimated introgression among lineages using additional genome data and found evidence for ancient hybridization leading to the common ancestor of Heliconius hortense and Heliconius clysonymus. We estimated the Ka/Ks ratio for each orthologous cluster and performed further tests to demonstrate genes showing evidence of adaptive protein evolution. Furthermore, we characterized patterns of expression for a subset of these positively selected orthologs using qRT-PCR. Taken together, we identified candidate eye genes that show signatures of adaptive molecular evolution and provide evidence of their expression divergence between species, tissues, and sexes. Our results demonstrate: 1) greater evolutionary changes in younger Heliconius lineages, that is, more positively selected genes in the cydno-melpomene-hecale group as opposed to the sara-hortense-erato group, and 2) suggest an ancient hybridization leading to speciation among Heliconius pupal-mating species.

Butterfly Mimicry Polymorphisms Highlight Phylogenetic Limits of Gene Reuse in the Evolution of Diverse Adaptations.

Some genes have repeatedly been found to control diverse adaptations in a wide variety of organisms. Such gene reuse reveals not only the diversity of phenotypes these unique genes control but also the composition of developmental gene networks and the genetic routes available to and taken by organisms during adaptation. However, the causes of gene reuse remain unclear. A small number of large-effect Mendelian loci control a huge diversity of mimetic butterfly wing color patterns, but reasons for their reuse are difficult to identify because the genetic basis of mimicry has primarily been studied in two systems with correlated factors: female-limited Batesian mimicry in Papilio swallowtails (Papilionidae) and non-sex-limited Müllerian mimicry in Heliconius longwings (Nymphalidae). Here, we break the correlation between phylogenetic relationship and sex-limited mimicry by identifying loci controlling female-limited mimicry polymorphism Hypolimnas misippus (Nymphalidae) and non-sex-limited mimicry polymorphism in Papilio clytia (Papilionidae). The Papilio clytia polymorphism is controlled by the genome region containing the gene cortex, the classic P supergene in Heliconius numata, and loci controlling color pattern variation across Lepidoptera. In contrast, female-limited mimicry polymorphism in Hypolimnas misippus is associated with a locus not previously implicated in color patterning. Thus, although many species repeatedly converged on cortex and its neighboring genes over 120 My of evolution of diverse color patterns, female-limited mimicry polymorphisms each evolved using a different gene. Our results support conclusions that gene reuse occurs mainly within ∼10 My and highlight the puzzling diversity of genes controlling seemingly complex female-limited mimicry polymorphisms.